Soft gluon emission off a heavy quark revisited 
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An improved generalized suppression factor for gluon emission off a heavy quark is derived within 
perturbative QCD, which is valid for the full range of rapidity of the radiated gluon and also has 
no restriction on the scaled mass of the quark with its energy. In the appropriate limit it correctly 
reproduces the usual dead cone factor in the forward rapidity region. On the other hand, this 
improved suppression factor becomes close to unity in the backward direction. This indicates a 
small suppression of gluon emission in the backward region, which should have an impact on the 
phenomenology of heavy quark energy loss in the hot and dense matter produced in ultra relativistic 
heavy ion collisions. 
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The main aim of ongoing ultra relativistic heavy ion 
collisions is to study the properties of nuclear or hadronic 
matter at extreme conditions. A particular goal lies in 
the identification of a new state of matter formed in 
such collisions, the quark-gluon plasma (QGP), where 
the quarks and gluons are deconfined from the nucleons 
and move freely over an extended space-time region. Var- 
ious measurements taken at CERN Super Proton Syn- 
chrotron (SPSHll and BNL Relativistic Heavy Ion Col- 
lider (RHIC) I^tSI lead to a wealth of information for 
the formation of the QGP during the first several fm/c 
of the collisions through the hadronic final states. New 
data from experiments at CERN Large Hadron Collider 
(LHC) have supported the existence of such a state of 
matter. 

Some of the important features of the plasma pro- 
duced in heavy ion collisions include energy loss and 
jet quenching of high energetic partons, viz., light and 
heavy quarks. The Gunion-Bertsch (GB) formula [l^l 
for gluon emission from the processes qq — > qqg has 
been widely used in different phenomenological studies 
of heavy ion collisions, in particular for radiative energy 
loss of high e nergy partons propagating through a ther- 



malized QGP |11| - |24] . The energy loss is presently a field 
of high interest in view of jet quenching of high energy 
partons, viz., bot h lig ht (22l . f25l - l27| and heavy quarks [lll - 



MllMa 12142a 128H321. Generally, one expects that jet 
quenching for heavy quarks should be weaker than that 
of light quarks. In contrast the non-photonic data at 
RHIC [sl reveal a similar suppression for heavy flavored 
hadrons compared to that of light hadrons. 

An early attempt to calculate the heavy quark energy 



loss in a QGP medium was done in Ref. [12| by using 
the GB formula of gluon emission for light quark scatter- 
ing [lo| and just modifying the relevant kinematics for 
heavy quarks. Later the soft gluon emission formula for 
heavy quarks in the high energy approximation [l3| was 



renewed in Ref. [13| for the small angle limit. Soft gluon 
emission from a heavy quark was found to be suppressed 
in the forward direction compared to that from a light 
quark due to the mass effect (dead cone effect). The 
corresponding suppression factor was obtained as [l^ , 

where Oq = M/E ^1. i? is the energy of the heavy 
quark with mass, M and 9, the angle between the heavy 
quark and the radiated gluon. Often in the literature 
I5L Til the expression 
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is used as the suppression factor to calculate heavy quark 
energy loss in heavy ion collisions. Here, k± denotes the 
transverse momentum of the emitted gluon and is related 
to its energy w by fcj^ = wsin^. For small angles (HI) re- 
duces to IT]). However, ^ produces not only a dead cone 
in forward direction (9^1) but also in the backward re- 
gion {9 ~ ±7r). The uniform use of such a dead cone 
for heavy quark energy loss may not be accurate enough 
since high energy scatterings off partons are associated 
with gluon emission in all directions with varying mag- 
nitude [T3l, as we will see below. 

In this article we revisit the issue and make an attempt 
to generalize the gluon emission off a heavy quark by re- 
laxing the constrains imposed in earlier calculations on 
the emission angle of the radiated gluons and the scaled 
mass of the heavy quark with its energy. We have found 
a generalized expression of the suppression factor that is 
identical to ([TJ and (I2|) for large E and small 9 but un- 
like ^ smoothly becomes unity (no suppression) in the 
backward direction. This supports the point of Ref. [13] 
that the main modification of the gluon radiation spec- 
trum due to a non-zero quark mass occurs at small angles 
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FIG. 1. Five tree level Feynman diagrams for the process 
Qq — >■ Qqg. In each diagram the thick upper line represents 
the heavy quark (Q) whereas the thin lower line represents 
the background light quark. 



(forward direction) and not at large angles (backward di- 
rection). 

In Fig. [1] the five Feynman diagrams for the process 
Qq — > Qqg are shown. According to the notation used in 
the figure, the Mandelstam variables are 



kA 



s = (fci + , s' = ih 

M = (fci - fc4)^ , u' = {k2 - K3 

t^{ki~ k^f , t' = (fc2 - 



with 



s + t + u + s' + t' + u'^ AA'P 



(3a) 
(3b) 
(3c) 



(4) 



Soft gluon emission (/cs — 0) ^33*] implies t' — >■ s' — )■ 
s, u' — u. In the center of momentum frame we consider 
the case where the energy of the emitted gluon, uj is much 
smaller than the momentum transfer yj\t\ « q±^ from the 
projectile (heavy quark) to the target (light quark) which 
again is small compared to the energy of heavy quark E. 
This leads to the hierarchy 



(5) 



It is important to note that the scaled mass of the heavy 
quark with its energy M / E and the gluon emission angle 
are free from any constrain. 

The gauge invariant amplitude for the process Qq ^• 
Qqg can be written as the squared matrix elements from 
the diagrams of Fig. [1] including their interference terms. 



'Qq9\ 
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where i and j run from to 4 and A^f^ = AiiAi* with 
Aii being the matrix element of diagram i (see Fig. [1]). 



With the hierarchy indicated in ([5]) the different matrix 
elements squared are obtained in the Feynman gauge as 
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For the calculation of the matrix elements we used RE- 
DUCE and CalcHep programs [34|]. We note that the 
results are in accordance with the matrix elements ob- 
tained from the process qq — ^ QQg [35'] by crossing the 
light antiquark and the heavy antiquark. Also, by proper 
transformation to light-cone variables, this result reduces 
to that obtained in Ref. within scalar QCD approx- 
imation and in light-cone gauge. 

The gauge invariant amplitude for the process Qq — > 
Qqg can now be obtained by summing all the sub- 
amplitudes ([7]), 



\MQq^Qgg\ ^ I2g'^ \MQq^Qq\ 



kl 



M2 ^ 



12g^\MQg^Q,\^^ll 



AP 



stan^d) 



^l2g^\MQ,^Q,\'^^[l + ^e'^^ 0) 

where 77 = — ln[tan(6'/2)], the rapidity of the emitted 
massless gluon. The two body amplitude is given by 



(10) 



Equation (|9]), which is the main result of the present 
article, carries a generalized suppression factor, 2? as 
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This improved suppression factor is valid in the full range 
of 9 (or rapidity of the emitted gluon) (i.e., —n < 9 < 
+Tr) and in the full range of M/^ (i.e., < M/^i < 1) 
as compared to Ref. [14]. As a note, the relation between 



3 



the center of mass energy y/s and the energy of the heavy 
quark E reads 



s = 2E^ + 2Ey/E^ - AP - 



(12) 



Below we discuss our results in more detail. First we 
consider two limits: 

1. Gunion-Bertsch limit: For M = Q, ^ reduces to 
the well known result of Gunion and Bertsch 
as 



\M 



qq'^qq'g\ — 12<? \Mqq^^qq'\ 
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where we have used (O that implies q± ^ k± [l9l - 



2. Dokshitzer and Kharzeev's result: In the limit M <Si 
and 6* < 1, it is ~ 2E and tan(6'/2) ~ 9/2 
and (O reduces to 



l^Qq^Qqq 



12g' \Mq,^q,\'^{1 



E^ 
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where 6*0 = M/E. This expression is precisely the 
result derived in Ref. [l3|. 

For convenience, we define TZ as the ratio of the squared 
matrix element of the 2 — >■ 3 to that of the 2 — > 2 pro- 
cesses, 



n = 



\M 



*Qq9\ 



„ 2 1 
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We note that this ratio, TZ is related to the gluon emis- 
sion multiplicity distribution 1^21| as dng/drjdk^ = 



TZ/IQ'K^. For the massless case, T^™-*'^ ig symmetric in 
rapidity. In contrast, a finite mass of the quark renders 
TZ to be asymmetric in rapidity. To explore this in more 
detail we consider the following rapidity regions: 

1 . Forward rapidity (?/ 3> 0) : In this case (|15p reduces 
to 



TZ. 



-2r, 



(16) 



Clearly, in this region of rapidity the gluon emis- 
sion is exponentially suppressed, which indicates 
the presence of the dead cone in the forward direc- 
tion if M ^ 0. 




FIG. 2. The suppression factor D from as a function of 
e and M/^. 



2. Mid-rapidity [r] Q): At mid-rapidity TZ depends 
only weakly on 77 as 



1 



AP 
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M2 



M2 



(17) 



3. Backward rapidity (77 <C 0): Here (jisp becomes 



(18) 



In this region the gluon emission does not depend 
on the mass and is, therefore, the same for heavy as 
well as light quarks. This is an important aspect for 
gluon emission off a heavy quark 12 , [l^ UsHitI [23| . 



We also note the dominant process (i.e., Qg — > Qgg) 
where a gluon acts as a target. Within the hierarchy ^ 
it differs from Qq — > Qqg only by a color Casimir factor 
Ca/Cf = 9/4 as 



since the two body part is given as 



(19) 



(20) 



and the other factors are the same for both processes in 
the considered approximations. Therefore, the factors T) 
and TZ remain unchanged. 

In Fig. [2] the suppression factor I) [cf. Eq. ([Tl]) ] is 
plotted as a function of 9 and M/ y/s. Around 6* ^ 1 we 
observe a canyon for small M/y/s and a valley for large 
Al/y/s, which clearly indicate a presence of a dead cone 
in the forward direction with respect to the propagating 
heavy quark. The spread of the dead cone increases as 
Al/y/s increases. In the backward region, 9 ~ ±7r, the 
suppression factor saturates to unity. This suggests that 
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FIG. 3. (color online) Comparison of our result Eq. (Illf) for 
the suppression factor I? to Eq. as a function of 6 for 
charm and bottom quarks with energy E = & GeV. As a 
note, a mass of 1.2 GeV (4.2 GeV) for charm (bottom) quarks 
has been used. 



the quark mass plays only a role in the forward direction 
when the energy of the quark becomes of the order of its 
mass. 

The possibility of this large angle scattering might be 
important for heavy-ion phenomenology in the context 
of the non- photonic electron data at RHIC and LHC. 
Furthermore, it might also have an impact on the de- 
scription of the forward-backward asymmetry of dijets 
and the seen ener gy d eposition at large angles in respect 
to the leading jet [36>| . 

Figure [3] compares our result for the generalized sup- 
pression factor T) in pT|) to that given in ([2]) as a function 
of the emission angle 6 for charm and bottom quarks with 
E = 6 GeV. Equation ([2]) agrees with our result in the 
domain of a small emission angle. However, the little 
variation in this region is due to the constrain M <C \/s 
employed in earlier calculations, whereas no such con- 
strain is set in our calculation. In contrast to ^ our 
result for the suppression factor V approaches to unity 
for large emission angles. This indicates that the back- 
ward emission is as strong as for light quarks. 

In summary, we derived a compact expression that con- 
tains a generalized suppression factor for gluon emission 
off a heavy quark through the scattering with a light par- 
ton. In the appropriate limit this expression reduces to 
the usually known dead cone factor. Our analysis shows 
that there is a suppression of soft gluon emission due to 
the mass of the heavy quark in the forward direction. On 
the other hand, the present findings also indicate that a 
heavy quark emits a soft gluon almost similar to that 
of a light quark in the backward rapidity region. This 
result might have important consequences for a better 



understanding of heavy flavor energy loss in heavy ion 
collisions. 
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